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Abstract 

The quasiclassical Green function formalism is used to describe charge and spin 
dynamics in the presence of spin-orbit coupling. We review the results obtained for 
the spin Hall effect on restricted geometries. The role of boundaries is discussed in 
the framework of spin diffusion equations. 
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1 Introduction 



A transverse (say along the y-axis) ^-polarized spin current flowing in response 
to an applied (say along the x-axis) electrical field 

Jl y = a sH E x , (1) 

is referred to as the spin Hall effect, and <j s h is called spin Hall conductivity [112] . 
This effect allows for the generation and control of spin currents by purely 
electrical means, which is a great advantage when operating electronic de- 
vices. Physically, the possibility of a non-vanishing a s n is due to the presence 
of spin-orbit coupling, which in semiconductors may be orders of magnitude 
larger than in vacuum. Clearly, two key issues are (i) how sensitive is the 
spin Hall conductivity to various solid state effects like disorder scattering, 
electron-electron interaction etc. and (ii) how the spin current can actually 
be detected experimentally. About the first issue, there is now a consensus 
that the effect of disorder scattering depends on the form of the spin-orbit 
coupling. In the case of the two-dimensional electron gas with a Rashba type 
of spin-orbit coupling, arbitrarily weak disorder leads to the vanishing of the 
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spin Hall conductivity in the bulk limit [3,4,5,6,7,8,9]. Concerning the second 
issue, the first experimental observations 1 1 |[TT] of the spin Hall effect have 
been achieved by measuring, optically, the spin polarization accumulated at 
the lateral edges of an electrically biased wire. Hence, the understanding of the 
spin Hall effect involves the description of boundaries. In order to address the 
two issues mentioned above, we develop in the following a quasiclassical Green 
function approach for the description of charge and spin degrees of freedom 
in the presence of spin-orbit coupling. 



2 The Eilenberger equation 



In this section we sketch how to derive the Eilenberger equation for the qua- 
siclassical Green function in the presence of spin-orbit interaction [12]. We 
consider the following Hamiltonian 

H=£ + b-*, (2) 

where b(p) is an effective internal magnetic field due to the spin-orbit coupling. 
For instance, in the Rashba case we have b = ap x e z . The Green function 
G tlt2 (xi, x 2 ) is a matrix both in Keldysh and spin space and obeys the Dyson 
equation (H — 1) 

id tl + —dl +fi- b(-i<9 Xl ) • <x) G tl t 2 {*i, x 2 ) = 5(h - t 2 )5(xi - x 2 ). (3) 
2m J 

A quasiclassical description is possible when the Green function depends on 
the center-of-mass coordinate, x = (xi + x 2 )/2, on a much larger scale than on 
the relative coordinate, r = xi — x 2 . In this situation, by going to the Wigner 
representation 

G(x 1 ,x 2 )=^e i P- r (5(p,x) (4) 
p 

and subtracting from Eq.Q its complex conjugate, one obtains a homogeneous 
equation for G(p,x) 

Sdt6+U^ + \d v {h-a),dJ^\-\h'tr,6\ = \Z,G], (5) 
2 2 J L j l j 

where only the center-of-mass time, t = (t\ + t 2 )/2, appears. On the right- 
hand-side of Eq.Q we have also introduced the self-energy 

EG(p,x), (6) 



2nN r - 

which takes into account disorder scattering at the level of the self-consistent 
Born approximation. In the spirit of the quasiclassical approximation, we make 
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the following ansatz for G 
G R G K 



G 



G J 




G R 



~ g R ~ g K 
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(7) 



where we assume that g does not depend on the modulus of p but at most 
on its direction. In this way we have separated the fast variation of the free 
Green function in the relative coordinate r from the slow variation of g in the 
center-of-mass coordinate x. The quasiclassical Green function is defined as 
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dC G(p,x), 



where £ = p 2 /2m — /i and p = p/p. Using the ansatz above, the ^-integration 
can be done explicitly. By assuming that b = 6(p)b(p) and 
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(9) 



(10) 



we find 

7T J-oo iV 

where N + and iV_ are the densities of state of the spin-split bands. In the 
absence of spin-orbit coupling, N± = N , and the function g coincides with 
the ^-integrated Green function g. In the present case, however, 



9 = ^o^}- 



By inverting Eq.( JTTT) one has 



9 = ^{(9o)~\9} + ^ 



1 r 



(9*a 



R\-l 



9o,9 



1 



N + -N_ 

2Nn 



b<r ),g , (12) 



where the last approximation, to be used in the following, holds for b <C ep- 
The ^-integration of G multiplied by a momentum dependent function leads 
to 

(13) 



1 



dtm(p)G=-{Mg£,g}. 

71 J-oo Z 

where M = m(p + )P + + m(p_)P_, and p± are the Fermi momenta in the two 
bands. To first order in b/ep one obtains 



^{Mg* g} « ^{M, g} = m(p + )g + + m(p_)g_ 



(14) 



where g± = §{P±, g}- Hence the standard procedure to obtain the Eilenberger 
equation for g via the ^-integration of Eq.flS]) yields 



E ( d t9u + 
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t,g . (15) 
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Eq. (|l5|) holds even for internal fields b for which the factorization b = b(p)b(p) 
is not possible, as long as one can assume |b| <C e F . For more details see |12j . 

From the above Eilenberger equation the expression for the charge and spin 
currents are readily obtained, since the Fermi surface average (...) of Eq.( fT5i) 
has the form of a continuity equation 



d t (g c )+dx-h = (16) 
d t (g x ) + <9 X ■ J: = 2 £ (b„ x g u ) x (17) 

v=± 

d t (g y ) + d x -jy = 2j2(Kxgv)y (18) 

v=± 

d t (g z ) + d*-J z s = 2j2(Kxgu) z , (19) 

v=± 



where we expanded the Green function in charge and spin components, g = 
g c + g • cr. For instance, the spin current with spin polarization along the e z 
axis reads 

AT f 

Zfrt) = --fJde[3fax,t)] K , (20) 
with J*(e;x,t) = (v F g z ). 



3 Spin Hall effect 



Before considering explicitly the role of the boundaries, it is useful to see how 
the spin Hall effect in the bulk may be analyzed in the present formalism. 
For a stationary and space independent case the equation for the Keldysh 
component can be written in the following way (a = 2ap F r) 
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In the above we have introduced the mean free path / = vpr. The presence of 
the electric field, along the e x axis, is accounted for by the source term Sg 




( 







S E = -4\e\Elf'(e) 
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+ — 



Pi-Pi 

V ) 



(22) 



which has been introduced in the Eilenberger equation by exploiting the gauge 
invariance in the derivative terms. In the above /(e) is the Fermi function. 
Notice that spin-charge coupling effects arise at the order of ajvp. From the 
angular average of Eq. (l2~T|) one observes immediately that (p x gf) = (P y g^) = 
0, i.e. no spin current with polarization along e z flows in the system. It is also 
instructive to first express gf in terms of the angle-averaged quantity (g K ) 
and then multiply by p y and take the angular average. The spin current reads 



When comparing to a calculation of the spin-Hall conductivity within the 
standard Kubo-formula approach, one realizes that the first term in square 
brackets corresponds to the simple bubble diagram, while the second term 
accounts for vertex corrections due to the spin- dependent part of the velocity 
operator. This second term is related to a voltage induced spin polarization 
in e y direction that was first obtained by Edelstein [15] ; s y = — \e\aTN E. 

In summary, we find that under the very special conditions we considered 
here the spin-Hall current is zero. We started from the Rashba Hamiltonian, 
that has a linear-in-momentum spin orbit coupling, we assumed that elastic 
scattering is spin-conserving, and finally we considered the stationary solution 
in the bulk of the two-dimensional electron system so that the spin density 
neither depends on space nor time. Relaxing one or more of these conditions 
finite spin-Hall currents are expected. 



4 Boundary conditions and diffusion equation 

The derivation of the boundary conditions for the quasiclassical Green func- 
tion is a delicate task, due to the fact that the boundary potential typically 
varies on a microscopic length scale which is shorter than the quasiclassical 
resolution. Therefore one must include the boundary effect in the very process 
of the derivation of the Eilenberger equation. This is especially important for 
interfaces between two different materials, where transport occurs through a 
tunnel junction or a point contact. Often a boundary with vacuum is simpler 



then 
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Fig. 1. Orientation of the boundary and the electric field. 



to describe since one couples only one incoming with one outgoing channel. 
In the presence of spin-orbit coupling, the non-conservation of spin and in 
particular the beam splitting, i.e., one incoming channel with direction pj n 
can be scattered into two outgoing channels p ou t, makes the problem more 
difficult and there is not yet a complete derivation of boundary conditions for 
the quasiclassical Green function. 

In the following we limit to two special types of boundary conditions and dis- 
cuss their role as far as spin relaxation and the spin Hall effect are concerned |14j. 
The first type of boundary conditions requires spin conservation, i.e. the spin 
currents normal to the interface are zero. We refer to them in the following 
as hard-wall boundary conditions, (see also [T5|rT6] ). For smooth confining po- 
tential it has been pointed out, that specular scattering involves some spin 
rotation in such a way to keep the scattered particle in the eigenstate of the 
incoming one [TTJ . We refer to this situation as soft- wall boundary conditions. 

The matching condition for the quasiclassical Green function can be obtained 
following the approach of Ref. [18] and reads 



where the matrix S describes the scattering at the boundary. 

To illustrate the role of boundaries, we consider the diffusive limit where a 
detailed analytical treatment is possible and which is also relevant for actual 
experiments. We assume that our two-dimensional electron gas is limited to 
the half-space y > and that there is a uniform electric field along the 
direction, see Fig. HJ Since the system is translational invariant along the 
direction, we consider only space dependence with respect to the e y direction. 
The diffusion equation reads 




(24) 



6 



d f -Dd: 



(d t - Dd 2 y ) 

(d t - Ddl) 



2Bd y s x 



1 

S x 

hs 

T s 

2 

s 

Ts 



y 



2Bd y p 

- s ) + 2Cd y s z 



2C8ySy, 



(25) 
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(27) 

(28) 



with D 



1 2 

2 V F r i r s 



r/[2(appT) ], B = Aa p F r , C = v F ap F r, and s 



— \e\ciTNoE is the bulk spin polarization in the presence of an electric field, 
mentioned at the end of the previous section. A general solution of the diffusion 
equations (I25H28I) can be found in the form 



s(y,t) 



(29) 



where a static solution requires 7 = 0. In the absence of electric field (i.e., 
so = 0), as in optical spin excitation experiments, the solution with longest 
lifetime has a finite wave vector in contrast to standard diffusion. The presence 
of boundaries does affect this result. With hard- wall boundary conditions 



- Dd y s x = n-j^O, (30) 
-DdySy - Cs z = n • j y = 0, (31) 
-Dd v s z + C{sy-s o ) = n-j g = 0, (32) 

where n is in the ^-direction, the mode with longest lifetime is localized at the 
boundary. With soft-wall boundary conditions, one has 

s x = and s y = so, (33) 

while the z-component of the spin is still conserved and therefore Eq. f|32l) 
remains valid. Due to the decoupling, at the boundary, of the three spin com- 
ponents, the mode with longest lifetime is no longer confined to the boundary 
and has s x oc sm(qy), s y — s oc sin(gy), s z oc cos(qy) with q ~ L" 1 (as in the 
bulk), where L s = \fDr s is the spin diffusion length. Although the two types 
of boundary conditions have a different effect in a time-dependent optical spin- 
relaxation experiment, both imply that, in the absence of an electrical field, 
the only static solution of the diffusion equation is with vanishing polarization 
(See Ref. [14J for further details). 

The presence of an electric field does not change this result provided the 
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polarization in the e y direction is replaced by the difference s y — sq, which 
enters both the diffusion equations (I25H28P and the boundary conditions (1301 - 
1551) . Hence, in the presence of an electric field, the only static solution has 
s y = s , s x = and s z = 0. Furthermore, the second condition in (T3HT) means 
that, with soft-wall boundary conditions and at the level of diffusive accuracy, 
the time-dependent approach to the static solution becomes infinitely fast at 
the boundary, in contrast to what happens in the bulk where relaxation occurs 
in a finite time. 



5 Numerical results 



In this section we give a brief overview of numerical results for both spin Hall 
effect and spin accumulation in finite systems. In the following we consider the 
Rashba Hamiltonian on a rectangular strip of length L x and width L y which 
is connected to reservoirs at x = and x = L x . We impose soft-wall boundary 
conditions at y — and y — L y . At t — the electrical field E is applied 
in x— direction. In order to solve the time-dependent Eilenberger equation 
(TT5l) numerically we introduce the discretization Ax and At for space and 
time, respectively. Our algorithm is exact to order (At) 2 and (Aa;) 4 and yields 
stable solutions for arbitrary long times provided At is chosen small enough. 
The boundary conditions are imposed such that at the border for incoming 
momenta the Green function is calculated from the Eilenberger equation while 
for outgoing momenta the scattering condition Eq. (T2~3| is applied. Fig. [2] shows 
the time evolution of the spin polarization on a line across the strip at x — 
L x /2, the Rashba parameter is a/vp = 10~ 3 . In the upper part of the figure 
the elastic scattering rate r is chosen such that appT = 1 which is already 
beyond the reach of the diffusion equation approach. The y— component s y of 
the spin polarization (left panel) increases monotonically from zero to the bulk 
value so while s z (right panel) shows strong oscillations close to the boundaries 
which decay after a few periods. These oscillations are associated with a spin 
current polarized in z— direction and flowing in y— direction. The magnitude 
of this spin current corresponds to a spin Hall conductivity <j s h — e/8ir which 
however persists only on a short time scale after switching on the electrical 
field. In the stationary limit only the bulk polarization of s y survives while 
the polarization of s x and s z is restricted to a small region around the corners 
of the strip [3]fT2] . In addition, the shape of these corner polarizations depends 
strongly on the details of the coupling between the reservoirs and the strip. 
In the lower part of Fig. [2] we choose appT = 0.1, i.e. the spin dynamics is 
diffusive. Note that at the boundaries s y approaches so on a much shorter time 
scale than in the bulk, in accordance with the boundary condition (1331) . Due 
to the overdamped dynamics of the spins there are neither in s y nor in s z time 
dependent oscillations. 
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Fig. 2. Voltage induced spin polarization as function of y and t at x = L x /2 on a 
strip of length L x = 201 and width L y = 101 for a/vp = 1CT 3 . The upper figures 
are obtained for apFT = 1, whereas in the lower figures apFT = 0.1. Left panel: s^, 
right panel: s z , both in units of the bulk polarization sq. 

To summarize, the quasiclassical approach provides a versatile theoretical 
framework for the description of the coherent spin dynamics in confined elec- 
tron systems in the presence of spin-orbit coupling. 

This work was supported by the Deutsche Forschungsgemeinschaft through 
Sonderforschungsbereich 484. 
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